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A THEORETICAL ANALYSIS OF THE EFFECT OF TIME LAG IN AN AUTOMATIC STABILIZATION 
SYSTEM ON THE LATERAL OSCILLATORY STABILITY OF AN AIRPLANE ' 

By Leonard Sternpield and Ordway B. Gates, Jr. 


SUMMARY 

A method is ‘presented Jor determining the effect of time lag 
in an automatic stabilization system on the lateral oscillatory 
stability of an airplane. The method is based on an analytical- 
graphical procedure. The critical time lag of the airplane- 
autopilot system is readily determined from the frequency- 
response analysis. 

The method is applied to a typical present-day airplane 
equipped with an automatic pilot sensitive to yawing accelera- 
tion and geared to the rudder so that rudder control is applied 
in proportion to the yawing acceleration. The results calculated 
tor this airplane-autopilot s^jstem by this method are compared 
with the airplane motions calculated by a step-by-step procedure. 

INTRODUCTION 

Recent calculations and. flight tests of several airplanes 
designed lor operation in the transonic speed range have 
indicated unsatisfactory damping of the lateral oscillation. 
The results presented in reference 1 show that the oscillatory 
stability can be improved by the use of an automatic pilot. 
The calculations of reference 1, however, were made on the 
assumption of an idealized control sj’stem without lag. 
Reference 2 points out that lag of the type in which the 
amount of control applied at a given instant is assumed to 
be proportional to a deviation which existed at a fixed time 
previous to the given instant can be represented mathe- 
matically by use of the lag operator where r, is 

nondimensional time lag based on the span and D>, is the 
differential operator. Lag of this type is generally referred 
to as time lag. Reference 2 suggests that for purposes of 
calculation the lag operator may be approximated by tkree 
terms of the Taylor’s series for T his approximation 

was used in some recent calculations (reference 3) and the 
results were found to be erroneous and misleading. The 
purpose of this report is to present a satisfactory method 
for determining the effect of time lag on the lateral oscilla- 
tory stability based on the exact expression of the lag 
operator rather than any approximation. Some recent 
analyses on the same problem, unknown to the authors at 
the time this problem was being analyzed, are presented in 
references 4 and 5. 

1 Supersedes NACA TN 2005, "A Thcoreticftl Analysis 
Sternfield and Ordway B. Gatos, Jr., 1950. 
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SYMBOLS AND COEFFICIENTS 
<t> angle of roll, radians 

^ angle of yaw, radians 

/3 angle of sideslip, radians {vjV) 

r yawing angular velocity, radians per second 

{d\}/jdt) 

yawing angular acceleration, radians per 
second per second {dll/flf) 

p rolling angular velocity, radians per second 

(d<j>ldt) 

0 sideslip velocity along lateral axis, feet per 

second 

V airspeed, feet per second 

P mass density of air, slugs per cubic foot 

q dynamic pressure, pounds per square foot 

(!'”■) 

b wing span, feet 

(S' wing area, square feet 

W weight of airplane, pounds 

m mass of airplane, slugs (Wjg) 

g acceleration due to gravity, feet per second 

per second 

Hh relative-density factor (mlpSb) 

y inclination of principal longitudinal axis of 

airplane with respect to flight path, positive 
when principal axis is above flight path at 
nose, degrees 

7 angle of flight path to horizontal axis, 

positive in a climb, degrees 

kx„ radius of gyration in roll about principal 

longitudinal axis, feet 

kz^ radius of gyration in yaw about principal 

vertical axis, feet 

Kx„ nondimensional radius of gyration in roll about 

principal longitudinal axis {kxjb) 

Kz^ nondimensional radius of gyration in yaw 

about principal vertical axis {kzjb) 

Kx nondimensional radius of gyration in roll 

about longitudinal stability axis 

{■yjKxf cos^ y-\-Kzf sin^ y) 
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K, 

Kxz 


C: 


Cn 


Cy 


^ d/3 

r 

dl3 


Cy = 


<7„ 


Cn, 




Cy = 


Cy, 


Cl- 


c.^r 

t 

Si 


drv 

. rh 
^2V 
(>C„ 

^PP 

2V 

dC, 

d PP- 

2V 

^Cy 

d Pp 

2V 

dCy 

^2V 

) P-P^ 

2V 

dC„ 

dSr 


Di 

P 

Ti/2 

Sr 

a 


Ws="p “ 

X = (I “H i coj 

T 


nondimensional radius of gyration in yaw 
about vertical stability axis 

{^JKzg^ cos^ v+Kx„^ sin^ rj) 
nondimensional product-of -inertia parameter 
((KzJ^—Kxg) sin V cos v) 

trim lift coefficient 

rr- ■ . /Rolling momentX 

rolling-moment coefficient ( ) 




J 


) 


ai ■ , /Yawing momentX 

yawing-moment coefficient^ ) 

, , ,, ^ /Lateral forceX 

lateral-iorce coefficient^^ ^ j 


time, seconds 

nondimensional time parameter based on span 
(Vtlb) 

differential operator 

period of oscillation, seconds 
time for amplitude of oscillation to damp to 
one-half its original value 
deflection of control, radians 
real part of complex root of characteristic 
stability equation 

angular frequency, radians per second 



time lag between signal for control and its 
actual motion, seconds 


(ts)c 

Kx 


V2 

Jl s ^ 2 ^ 


critical time lag 

maximum amplitude of acceleration in yaw 
produced by control deflection of unit 


amplitude ( — 


(I ) 




amplitude of control-surface oscillation pro- 
duced by autopilot in response to oscillation 

1 5. 


of airplane acceleration 


(-- ) 

\ Ip autniiiloi r 



(Kc)s — 


ks 


e 

6 a 


phase angle, radians 

phase angle of lag of 5^ behind ip when oscillat- 
ing control surface forces airplane to 
oscillate, radians 

phase angle obtained from frequency response 
of autopilot 


ANALYSIS 

The investigation of the effect of time lag on the lateral 
oscillation may be conveniently divided into two parts: 
(a) determination of the smallest time lag which would result 
in a neutrally stable oscillation, referred to as critical time lag, 
and (b) the effect of a given time lag on the lateral oscillatory 
stability. It is important to know some of the results of the 
analysis obtained in part (a) in order to facilitate the analysis 
presented in part (b). Part (a) is based on the frequency- 
response method of analysis (references 6 and 7). This 
method affords a relatively simple means of determining the 
critical time lag of an automatic stabilization system and 
thereby of establishing the range of time lags for which the 
airplane motion is stable. Part (b) treats the solution of a 
transcendental equation by means of an analytical-graphical 
procedure. "The analysis is presented for an airplane 
equipped with an automatic pilot sensitive to yawing ac- 
celeration and geared to the rudder so that rudder control 
is applied in proportion to the yawing acceleration, as sug- 
gested in reference 3. A similar analysis is applicable, how- 
ever, to any automatic stabilization system with time lag. 

The equations of motion used are expressed in terms of 
the nondimensional time parameter based on the span of the 

airplane Si=^> but the results of the calculations obtained 

in terms of Si have been converted from the nondimensional 
time Si to time t in seconds. Thus the discussion of the 
results and the figures included in the report are given in 
terms of t. 

DETERMINATION OF CRITICAL TIME LAG 

The critical time lag of a system is defined as the time lag 
that results in a neutrally stable or steady-state oscillation. 
The motion of the control 5, and the airplane acceleration \p 
when a critical time lag exists are shown in figure 1. This 
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Figure 1.— The relationship between the airplane acceleration and the control motion for the 


case of critical time lag. 0 =wA?; Ai=T=— • 


figure indicates that the relationship between the time lag t 
and the phase angle 6 between the motion of the control and 

Q 

the airplane acceleration can be expressed as t= — Since 

CO 

r can be expressed as a function of 6, the frequency-response 


method can be utilized to determine the critical time lag 
for a combined ahplane-autopilot system. This method re- 
quu’es that the frequency-response curves be calculated for 
the airplane and autopilot separately and the results analyzed 
to ascertain the conditions for neutral oscillatory stability 
of the airplane-autopilot system. 

Frequency-response curves for airplane. — The frequency- 
response curves for the airplane are obtained from the cal- 
culation of the steady-state motion of the ah'plane in response 
to a sinusoidal forcing function of unit amplitude (see ref- 
erence 8). Thus if 

6 r = sin Us-Sj, 

the acceleration of the airplane is 

D^^={Ka)s sin (cojSs-f 0^) 

The values of {Ka)s, known as the amplitude ratio, and 6a 
are obtained over the desired range of angular frequencies by 
substituting iojj for in the expression for which is 

derived from the lateral equations of motion. The non- 
dimensional lateral equations of motion, referred to stabilitj* 
axes, for a given control deflection arc 


^ 'I ' — ’■ (1) 

^ — 9 Ih — f — 2 (^y\Di , — (\ tan y^\l/-\-(2nt,Dti d = 0 


The dei ivativos fV 

®r 


d6r 


and 



dfv 


are usually very small and therefore have been neglected in equations (1). 


Hence, 






‘^lybKxzI^b" ^ 0„^Db 


0 


-a 


-\Cy^Db-CL 0 


5r 

2fyj<xWb^-^r\jA 

2f.bKxzDb^-^ CtDb 

-A 






-^Cy^Db-C^ 

^2/ii,— i Cy^ Db—CL tan 7 



( 2 ) 


After the numerator and denominator are expanded by tbe 
method of determinants, the expression for results in 

the ratio of two polynomials in The substitution of 

ico, for in equation (2) gives a complex number A+iB 
which may be expressed as (Ka)s6‘''^. The amplitude ratio 


5 , which is equal to Ka^> can be determined from the 

relation (Ka)s= -\] A?-\-B^. The phase angle 6a can be de- 

B 

termined from the relation 0 A = tan“^-r- 
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Frequency response of autopilot. — The frequency response 
of the autopilot is obtained from the equation for the control 
motion with time lag taken into account 

Sr = ksD,?il/{8„—Ts) (3a) 


where the term Dt^i}/{Sb — Tj) signifies the fact that the amount 
of control applied at a given instant is proportional to the 
acceleration at a fixed time Ts previous to the given instant. 
This time-lag effect can be expressed by the so-called lag 
operator Thus equation (3a) becomes 

(3b) 


Solving equation (3b) for gives 


where 


P.V 1 

Sr ks 


e's'^b={Kc)se"^°>> 


(TT \ ^ TV ^ 


autopilot 


(4) 


Substituting iug for Dt in equation (4) results in an amplitude 
ratio (Kc)s=^’ which is independent of frequencj', and a 

rC g 


phase angle dc=TsO>s- 

Conditions for stability. — The necessary and sufficient 
conditions for any one of the oscillatory modes to be neutrallj' 
stable are that, at a particular frequency, the phase angles 
and amplitude ratios of the airplane and autopilot must be 
equal — that is, 6 a — 6c and {Ka)!,={Kc)s- Even though 
these conditions are satisfied, the resultant motion of the 
airplane which is composed of all the individual modes of 
motion may be unstable since, as is pointed out subsequently, 
an additional unstable oscillatory mode may be present. In 
general, the condition which must be satisfied in order that 
all the oscillatory modes be stable is that, at each angular 


frequency where Oa = 6c, the ratio 


(Ka). 


must be less than 1 


VQ. 111^ ICClTlW /Tjf' \ 

at that frequency. The mathematical proof of this state- 
ment is given in reference 9. 

Illustrative example. — The foregoing method is applied 
to a typical present-day high-speed airplane having the 
characteristics presented in table I. The value of the control- 
gearing ratio k is arbitraril}^ assumed to be 0.0427. (This 
value of A:=0.0427 corresponds to a rudder deflection of 1° 
for a yawing acceleration of 23.4°/sec/sec.) The amplitude- 
ratio and phase-angle curves for the airplane, plotted as a 

V 

function of angular frequency co, where co=<»!s -j-i are pre- 


sented as solid-line curves in figures 2 (a) and 2 (b), respec- 
tively. As the frequency increases to infinity, Ka approaches 
a value of 15.98 and 6a approaches x. The dashed line in 
figure 2 (a), which is independent of frequency, is the ampli- 
tude ratio of the autopilot Kc. The phase-angle curves of 
the autopUot are straight lines with slopes equal to r, where 


T=y Ts, and are shown as dashed lines in figure 2 (b) for 

several values of t. Since the phase angle remains between 
the range of 0 to 2x, 6c=tw continues to repeat itself when- 
ever tco^2x. To take account of this fact, 6 c is plotted as 
a series of parallel lines for each value of t. Figure 2 (a) 
indicates that Ka = Kc at a>=3.8 and co=8.5. The corre- 
sponding values of r where 6a — 6c at a; = 3.8 and co=8.5 are 
t=1,63 and t=0.38, respectively. One of the oscillatory 
modes of motion is thus neutrally stable when t=1.63 and 
t=0.38. However, as mentioned previously, the motion of 
the airplane is neutrally stable only if all other oscillatory 
modes present are stable. This condition is satisfied for 
t=0.38, since for each value of co where 6a = 6c the ratio 

^<1. When r=1.63, one of the oscillatory modes is 

neutrally stable but the system is unstable, because at co=6, 

6a = 6c but ^^>1. An analysis indicates that for values 

of iircr<Cl5.98, which is the limiting value of Ka, the system 
will be unstable for any infinitesimal time lag. The reason for 
the instability is that, for an}^ infinitesimal time lag, the 
value of 6a is equal to 6c at some high frequency where it 

can be shown that^>l since at the very high frequencies 

Ac 

Ka {Ka)^^ + ^Ka 
K c (_Ka)w->«> KKc 

where AKa and AKc arc small incremental values. 

TABLE I.— STABILITY DERIVATIVES AND MASS CHAR- 
ACTERISTICS OF A TYPICAL PRESENT-DAY AIRPLANE 

W/S, lb/ft2 

5, fU 

6, ft 

p, slugs/fO 

V, ft/sec 

7 , deg 

Cl 

Ml. 

AV 

Kz^ 

Kxz 

1 ), deg 

Ci^, per radian. 

Cl , per radian. 

C„^, per radian 
per radian. 

Cy^ , per radian 
Cr^, per radian 
Cy^, per radian. 

per radian. 

Cip, per radian. 

C„j , per radian 


65 

130 

28 

0.00089 

797 

0 

0.23 

80.7 

0.00967 

0.0513 

-0.00145 

- 2.0 

-0.40 

0.08 

-0.0155 

-0.40 

0 

0 

- 1.0 

0.25 

-0.126 

-0.163 


ft radians 


EFFECT OF TIME LAG IN STABILIZATION SYSTEM ON LATERAL OSCILLATORY STABILITY OF AIRPLANES 


5 



(a) Amplitude ratio. 

(b) Fbase angle. 

Figure 2.— Frequency-response curves of the typical presont-day airplane and the assumed autopilot. 
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In order to verify the results predicted by the preceding 
analysis, motions of the airplane were calculated subsequent 
to an initial disturbance of 5° in sideslip. The calculations 
involved a step-by-step procedure based on the Kutta 3/8 
method (reference 10). The results of these calculations for 
t=0.38 and t=1.63 are presented in figures 3 and 4, respec- 
tively. The solid-line curve in figure 3 . was obtained by 
using a time increment of 0.095 and the motion is seen to be 
slightly unstable, whereas neutral stability is predicted by 
frequency-response analysis for t=0.38. An additional cal- 
culation was made by using a time increment of 0.0475, 
represented by the dashed-line curve in figure 3, and although 
the motion was still slightly unstable, the trend indicated 
by reducing the time increment was such as to make the 
oscillation more nearly neutrallj^ stable. The airplane mo- 
tion for T= 1.63 is presented in figure 4 and, as was predicted, 
the motion is unstable. A neutrally stable oscillation was 
also predicted for this value of time lag but it is apparent 


from figure 4 that the unstable mode influences the airplane 
motion more than the neutrally stable mode. 

EFFECT OF TIME LAG ON LATERAL OSCILLATORY STABILITY 

Derivation of equations. — The nondimensional equations 
of motion, referred to the stability axes, which include the 
effect of an autopilot applying rudder control in proportion 
to the yawing acceleration at time Sj— r„ are obtained by 
combining equation (3b) with equations (1). When 
is substituted for <f>, for xf/, and for /3 in the re- 

sultant equation written in determinant form, X must be a 
root of the characteristic stability equation 

.^X^ -b + CX^ -b Z?X + £’+ -f (7 V -f D'X) = 0 

(5) 

where A, B, C, D, E, A' , B' , C", and D' are functions of the 
mass and aerodynamic parameters of the airplane. The 



Kicvhe 3.--KfToct of r=0.38 on the motion of the airplune. 



EFFECT OF TIME LAG IN A STABILIZATION SYSTEM ON LATBBAL OSCILLATORY STABILITY OF AIRPLANES 7 


expi’essions for A, B, C, D, and E&ra given in reference 1 and 


B' — (2/j.i,Kx^O j 


C'=(-ic,C,,+ic,,0,,)c.._ 

Z?'=C7iC,,<7„,^ 


The damping and period of the lateral oscillation in 
seconds are given by the expressions 


Ty, 


-0.693 

V 


a 


p_^ A 

0 ,, V 


( 6 ) 


where a and co* are the real and imaginary parts of a complex- 
root of equation (5) . 

Determination of roots of transcendental equation. — The 

characteristic stability equation of this system (equation 
(5)) is seen to be a transcendental equation because a con- 
stant time lag in the automatic stabilization system is rep- 
resented by the so-called lag operator It is apparent 

that the complex roots of such an equation cannot be deter- 
mined by conventional methods. A method of obtaining 
the complex roots of this transcendental equation for a par- 
ticular value of Ts is therefore presented. 

If a-biwj is substituted for X in equation (5) and the real 
and imaginary quantities arc separated, two equations in a 
and coj result: 


■*“ = /*’, (a, cos) 

(7a) 


(7b) 




8 


REPORT 1018 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


A simultaneous solution of equations (7a) and (7b) gives the 
desired values of a and which satisfy equation (5). The 
method used to solve equations (7a) and (7b) simultaneously 
is basically a graphical procedure. For a series of values of 
o!s, the right-hand sides of equations (7a) and (7b) are 
plotted against a as illustrated in figure 5. The solid lines 
and dashed lines correspond to the functions Fi(a,co,) and 
F 2 (a,o>s), respectively. The variation of with a is also 
plotted in figure 5. The exact values of a and co, for which 
Fi(a,Ws)=F 2 (a,cos)==«“^"“ are determined from a cross plot 
of the results of figure 5 as shown in figure 6. In this figure, 
with a as the abscissa and Wj as the ordinate, are plotted the 
values of a and oj, which correspond to the intersection of 
the curve with the functions Fi(a,Ug) and i^ 2 (a,“s)- 

The solid curve in figure 6 thus satisfies the equation 
e-T,a—Fi(a,ois) and the dashed curve satisfies the equation 
^-r,a—F^{a,ois). The values of a and ojj at the intersection 
of these two curves therefore determine a root of the char- 
acteristic stability equation (equation (5)). The period and 
damping of the lateral oscillation are determined from equa- 
tions (6) by using this root. 



I I I I I I I I 

0 -a 


Figure 5. — The effect of a and on the functions Fi(a,w*), F 2 (a,wt), and 
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Figure 6.— A cross plot of the points of intersection of and f' 2 (a,w,) with e-%<»as 

determined from figure 5. 


An alternative method for the simultaneous solution of 

^ 2^2 ^ 3fj3 

equations (7a) and (7b) is to substitute 1 — Xsa-l — ^ 

for in each of these equations. The calculations in- 

volved become considerably less laborious when this ap- 
proximation is made and results obtained by using this alter- 
native method have been found to be in excellent agreement 
with results obtained from the method described previously. 
Figure 7 shows that for |r.,a| = l very close agreement is 
obtained between the exact value of and its approxi- 
mation by three and four terms of the series for In 

actual practice the product r^a will almost invariably be 
much less than 1. When this substitution is made, both 
equations (7a) and (7b) become seventh degree in a for a 
given value of cos. Although these equations are of high 
degree, no serious problem is presented since the values of 
a desired must be real and, in general, small; that is, only 
one or two of the roots of these high-degree equations are 
of interest. In order to calculate the complex roots of 
equation (5) for a particular value of the following pro- 
cedure should be used. For a sequence of values of Wj, 
compute a from equations (7a) and (7b). The results ob- 
tained from each equation may then be plotted in a figure 
similar to figure 6 and the complex root of the characteristic 
stability equation determined from the intersection of the 
two resulting curves. If the value of the approximate 
series for is in good agreement with the exact value 
of «“’■>“ for the a determined from the intersection of the 
two curves, then the complex root obtained is valid. A 
point of intersection for a value of a which would not give 
satisfactory agreement between the approximate series and 
might exist, however. If such be the case, the accurate 
point of intersection is readily ascertained. The estimated 
values of a and oig represent the point of intersection of 
Fi(a,ois) or F 2 {a,ois) with the series approximating 
The desired point is the point of intersection of Fi{a,Ws) or 
F 2 {a,o}s) with The first step is to evaluate the term 

for several values of a in the vicinity of the estimated 
point of intersection. The expressions F’i(a,a>s) and F 2 (a,Ws) 
are then evaluated for values of co^ slightly less than and 
greater than the estimated value of coj for several values of 
a in the range of the estimated point of intersection. Thus, 
the corrected curves of equations (7a) and (7b) are obtained 
and at their point of intersection, the accurate values of a 
and CO, are determined. 

A method for constructing curves of constant period 
and damping as a function of tj and ks is presented in the 
appendix. 

Range of co* and a to be used in the determination of the 
complex roots of the characteristic stability equation. — In 
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Ficurk 7. — A comparison between and Us approximation by threo and four terms of the series for 


general, the analysis required to determine the damping of 
the airplane-autopilot system would be carried out for those 
values of tj that result in a stable sj^stem, since the purpose of 
equipping an airplane with an autopilot is to stabilize or 
increase the damping of the airplane motions. The range 
of Ts for which the system is stable is determined from the 
frequency-response curves of the airplane and autopilot as 
indicated in the section entitled “Determination of Critical 
Time Lag.” Thus, if the system is stable, only negative 
values of a need be investigated. Also, if two stable oscil- 
latory modes of motion exist, the least stable one is of 
greatest interest — that is, the complex root with the smallest 
real part is the one of most importance. 

The estimated range of values for that should be used 


in the analysis to determine the effect of some particular 
value of time lag, located between Tj= 0 and the critical 
time lag (tj)c, on the damping of the oscillation is obtained 
from the known values of the frequency of the oscillation for 
the cases where Tj= 0 and Ts=(tj)c. The imaginary part of 
the comple.x root of equation (5), which becomes a quartic 
equation when 7-j=0, gives the value of coj for the case of 
Ts=0. The value of coj for Ts=(ts)c is determined from the 
analysis presented in the section entitled “Determination of 
Critical Time Lag.” As mentioned in reference 11, the 
frequency of the oscillation decreases as time lag increases; 
thus for a value of 0 <[ts<^(ts)c, the estimated values of coj 
should include frequencies greater than the value of co, at 
'Ts=(ts)c and less than the value of oj* at Ti=0. 
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(■d) r = 0. (<I) r=0.25. 

(b) T=0.10. {(?) t=0.287. 

(c) t=0.2U. (f) t=0.38 (same as fig. 3). 

FiouifE 8. — Effect of time lag on the airplane motion in sideslip. 
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Airplane motions for various values of t. — For puiposes of 
comparison, the motion of the airplane in sideslip subsequent 
to an initial displacement in sideslip of 5° was calculated, with 
the use of a step-by-step procedure, for t=0, 0.10, 0.20, 0.25, 
0.287, and 0.38. The results are presented in figures 8 (a) 
to 8 (f). These figures indicate that as r increases from 0 to 
0.2 the period increases slightly, whereas the damping is 
markedly improved. The frequency of the oscillation is 



(a) Damping. • 

(b) Period. 

Figurk 9.— Effect of time lag on the period and damping of tho lateral oscillation. 


about 01=3.7. However, as t continues to increase, the 
presence of a high-frequency oscillation is noted in the motion 
since this oscillation becomes less damped and the low- 
frequency oscillation becomes more heavily damped. (See 
figs. 8 (d) and 8 (e).) For r=0.38, the motion is neutrally 
stable at the high frequency of to =8. 5 but the low-frequency 
oscillation does not appear in this motion since it is very 
well damped. In figures 9 (a) and 9 (b), the period and 
damping of the lateral osciliation for several values of t, 
calculated by the method discussed in this report, are com- 
pared with the period and damping readily obtained from the 
motion calculations shown in figures 8 (a) to 8 (f) . The very 
good agreement between the results presented in figures 
9 (a) and 9 (b) would probably be improved if the increment 
selected for the step-by-step calculations were reduced. The 
trends indicated by these results, however, are applicable 
only to the particular airplane-autopilot system considered 
in this example and may not be generalized to any other 
arbitrary airplane-autopilot system. 

CONCLUDING REMARKS 

.\ method is presented for determining the effect of time 
lag in an automatic stabilization s 3 "stem on the lateral 
oscillator 3 ' stability of an airplane. The method is applied 
to a typical present-da 3 ^ airplane equipped with an automatic 
pUot sensitive to yawing acceleration and geared to the 
rudder so that rudder control is applied in proportion to the 
yawing acceleration. The results calculated for an airplane- 
autopilot system by the method described are in good 
agreement with the airplane motions calculated by a 
step-b 3 '-step procedure. 


IjAnglky Akron.\utical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 28, 1949. 


APPENDIX 

A METHOD FOR CONSTRUCTING CURVES OF CONSTANT PERIOD AND DAMPING AS A FUNCTION OF t. AND k, 


If an over-all picture is desired of the effect of time lag 
for a given value of ks or the effect of varying k, for a given 
time lag, construction of curves of constant period and damp- 
ing with the use of the following method is recommended. 

The characteristic stability equation of the aii’plane- 
autopilot system (equation (5)) may be rewritten in the 
form: 

If X=a-|-fc«)s is substituted in each side of equation (Al), 
the condition that X be a root of the characteristic equation 
is that the complex number Ai-\-iBi obtained from the left- 
hand side must be equal to the one obtained from the right- 
hand side A.2-t-'i.S2- The quantities Ai-\-iBi and 
may be represented by the expressions and 

respectively. Therefore, this requirement is equivalent to 
saying that and if ^ is to bo a root of equation 

(Al). If X=a-j-iws is substituted in the right-hand side, 
the following expressions result: 

^’se~’'*“(cos TsCOs — i sin Ts1^s)= A2AiBo=R2e’^2 

where 

and 

n 1 TsOis o 

C72=tan = — = — TsO>s 

cos TsO^s 

Therefore, if X=a-f-'i^s is substituted in the left-hand side and 
Ai~\-iBi=Rie'^i is obtained, the value of r required to make 
61 = 62 can be determined. Since Tj is therefore determined 
and a is fixed, the value of k^ necessary to make = can 
be calculated. Thus for these values of ks and t^, 

X = a f Ws 

isarootof the transcendental stability equation (equation (5)). 

For a given value of a, the analysis may be made through- 
out the range of and the corresponding values of ks and 
Ts determined. A curve representing this value of a may then 
be plotted as a function of ks and t,. This procedure is 
repeated for a sequence of values of a and the corresponding 
curves in the ks,Ts plane are plotted. Each point on a curve 
of constant a represents a particular value of cos. Curves of 
constant frequency may therefore be plotted by drawing a 
curve through the given value of coj on each one of the a 
curves. The values of a and Wj are converted to and P 
by equations (6) . The final result would consist of curves of 
constant damping and curves of constant period P in the 
ks,Ts plane. Thus the effect of time lag on the lateral oscil- 
lations for any value of the gearing ratio ks, or the effect of 
varying ks for any value of time lag, may be ascertained. 
For purposes of illustration several lines of constant Th and 
P in the k,r plane were calculated for the typical present- 
day airplane described in table I and are shown in figure 10. 

A curve can be plotted in the ks,Ts plane which divides the 
quadrant into a satisfactory and an unsatisfactory region 
12 
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Figure 10.— Several curves of constant period and constant damping for the typical presenf- 
day airplane described in table I. 

according to any prescribed relationship between the period 
and damping of the lateral oscillation. In order to calculate 
this curve, several values of a and Us that exactly satisfy the 
criterion should be selected and substituted in equation (Al). 
The combination of ks and ts is then obtained for each set of 
values of a and and the desired curve is plotted in the 
ks,Ts plane. 
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